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ABSTRACT 

While the importance of dusty asymptotic giant branch (AGB) stars to galactic chem- 
ical enrichment is widely recognised, a sophisticated understanding of the dust forma- 
tion and wind-driving mechanisms has proven elusive due in part to the difficulty in 
spatially-resolving the dust formation regions themselves. We have observed twenty 
dust-enshrouded AGB stars as part of the Keck Aperture Masking Experiment, re- 
solving all of them in multiple near-infrared bands between 1.5/im and 3.1/im. We 
find 45% of the targets to show measurable elongations that, when correcting for the 
greater distances of the targets, would correspond to significantly asymmetric dust 
shells on par with the well-known cases of IRC +10216 or CIT 6. Using radiative 
transfer models, we find the sublimation temperature of T su b (silicates) = 1130 ± 90 K 
and T su b (amorphous carbon) = 1170 + 60 K, both somewhat lower than expected from 
laboratory measurements and vastly below temperatures inferred from the inner edge 
of YSO disks. The fact that O-rich and C-rich dust types showed the same sublimation 
temperature was surprising as well. For the most optically-thick shells (r 2 . 2 gm > 2), 
the temperature profile of the inner dust shell is observed to change substantially, an 
effect we suggest could arise when individual dust clumps become optically-thick at 
the highest mass- loss rates. 

Key words: radiative transfer — instrumentation: interferometers — circumstellar 
matter — stars: AGB stage — stars: dust shells 


1 INTRODUCTION 

One of the most dramatic phases in the life of an inter- 
mediate mass star is the Asymptotic Giant Branch, a rel- 
atively short period where a star loses most of its initial 
mass through a dusty wind. Researchers still do not un- 
derstand all the ingredients necessary for producing the 
high mass- loss rates observed during this stage. The mas- 
sive envelopes ejected during this phase are thought to be 
later illuminated during the planetary nebula stage, a stage 
where most stars show strong bipolar circumstellar struc- 
tures (Balick & Frank 2002). 

Following the advent of infrared detectors, early workers 
made simple spherically-symmetric models of dusty shells 
around large samples of AGB stars fitting only to the spec- 
tral energy distributions^. g., Rowan- Robinson & Harris 
1982, 1983b, a). M-type stars are typically surrounded by 
dust shells composed of amorphous silicates while C-stars 
have carbonaceous dust. These early workers were able to 


show that dust condensed around 1000 K within a few stel- 
lar radii of the stars and also estimated mass-loss rates typi- 
cally 10 _6 M©/yr and as high as 10 -4 M©/yr. More recently, 
Ivezic & Elitzur (1995) developed the code DUSTY to study 
dust shells in a systematic way and made models for a large 
sample of stars, again fitting just the spectral energy distri- 
bution. 

The simple picture of spherically-symmetric and uni- 
form mass-loss was challenged by the observations of the 
Infrared Spatial Interferometer (ISI), a long-baseline mid- 
infrared interferometer (Danchi et al. 1994). These workers 
found a diversity of shell morphologies with some red gi- 
ants showing episodic dust shells ejections and others with 
a more continuous distribution of dust. A more dynamic 
and asymmetric vision of mass-loss fit into debates into the 
origins of bipolar symmetry in planetary nebulae. High an- 
gular resolution near-infrared speckle and aperture mask- 
ing on 8-m class telescopes were able to image fine de- 
tails on some dust shells, such as the prototype carbon 
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star IRC +10216 (e.g. Tuthill et al. 2000a). An elaborate 
model was presented by Men’shchikov et al. (2002) arguing 
for complex, spatially- varying dust properties and density 
structures. While IRC +10216 shows complexity within the 
inner few stellar radii, it is unclear if these structures rep- 
resent global asymmetries or just weather conditions of the 
dust formation process observed in situ. 

Here we present the full dataset of dust-enshrouded 
giants observed with the 10-year project called the Keck 
Aperture Masking Experiment (Tuthill et al. 2000b). This 
experiment delivered well-calibrated spatial information on 
the scale of ~50 milliarcseconds (mas) in the astronomical 
K band (Ao = 2.2/xm), enough to resolve all the dusty tar- 
gets presented here and to measure their dust shell sizes and 
asymmetries. This paper includes 20 objects with observa- 
tions in typically 3 wavelengths ranges, 1.65/im, 2.2/xm, and 
3.1/xm. We have also extracted photometry to construct co- 
eval near-IR spectral energy distributions - an important 
factor since these objects pulsate and show large variations 
in flux on yearly timescales. Lastly, we used a radiative 
transfer code to fit each epoch of each target star using 
simultaneously the NIR photometry and multi-wavelength 
angular size information from Keck masking. 

The primary goals of these observations and modeling 
efforts are to measure the physical characteristics of a large 
sample of the most extreme dusty AGB stars, to address the 
question of the onset of circumstellar asymmetries, to deter- 
mine any differences between silicate and carbon-rich dust 
shells, and to constrain the optical properties of the dust 
particles themselves. Lastly, this publication marks the final 
large data release of AGB star data from our diffraction- 
limited Keck masking experiment and we anticipate this 
work will provide a rich dataset for more detailed modelling 
efforts by other workers. 

2 OBSERVATIONS 

2.1 Overview of Observations 

Our observations consist of photometric and visibility data 
taken on 20 different stars at the W.M. Keck observatory 
between December 1997 and July 2002. The wavelengths 
at which these stars were observed and the properties of 
the corresponding Liters are listed in Table 1. A listing of 
the observed stars, segregated into carbon-rich and oxygen- 
rich groups, along with their basic properties can be found 
in Table 2. Most stars were measured at more than one 
epoch during this time span allowing for robust internal data 
quality checks. 

2.2 Photometric Data 

Aperture masking procedures consist of alternating target 
and calibrator observations that allow for basic photome- 
try in most observing conditions. As part of the standard 
pipeline (Monnier 1999; Tuthill et al. 2000b) we performed 
aperture photometry on each object, allowing the difference 
in magnitude (Amag) between the target star and calibra- 
tor star to be measured. The Vizier catalog service, most 
often referencing the Catalogue of Infrared Observations 
(Gezari et al. 1999) and 2MASS (Cutri 2003), was used to 


determine magnitudes at infrared wavelengths for the cali- 
brators. Interpolation was used between wavelengths found 
in the catalogues and the wavelengths at which our data was 
taken. Occasionally no mid-IR measurements were available 
for some calibrators and we used the calibrator spectral type 
and the K band flux to estimate the flux density at these 
longer wavelengths. 

As a data quality check we compared our photometry 
with 2MASS and found good general agreement, although 
strict agreement was not expected since our targets are 
highly variable and there is some difference in beam sizes. 
We estimated the error on the photometry points at 10% 
based on night-to-night variations. However, there were in- 
stances when we assigned larger errors (between 10 and 32%) 
due to saturation of the 2MASS photometry used for the 
calibrator, intrinsic variability of the calibrator, or effects 
of cirrus clouds in some of the original data. Indeed, there 
were some nights too contaminated by variable clouds to 
allow photometry to be extracted at all. 

Table 3 is a journal of observations, including the ob- 
serving date(s), the filter (s) used, the aperture mask(s) used, 
and calibrator star name. We have compiled the adopted 
calibrator properties in Table 4. 


2.3 Visibility Data 

2. 3. 1 Methodology 

Our group carried out aperture masking interferometry at 
the Keck-1 telescope from 1996 - 2005. We have published 
images and size measurements with (at the time) unprece- 
dented angular resolution on topics ranging from young 
stellar objects, carbon stars, red supergiants, and photo- 
spheric diameters of Mira variables (e.g., Monnier et al. 
1999; Tuthill et al. 2000a, b; Danchi et al. 2001). 

The NIRC camera with the image magnifier 
(Matthews et al. 1996) was used in conjunction with 
the aperture masking hardware to create fringes at the 
image plane. The data frames were taken in speckle mode 
(Tmt=0.14s) to freeze the atmosphere. In the work pre- 
sented here, multiple aperture masks and bandpass Liters 
were employed. After Lat-Lelding, bad pixel correction, and 
sky-subtraction, Fourier methods were used to extract fringe 
visibilities and closure phases from each frame and averaged 
in groups of 100 frames. Absolute calibration to account 
for the optical transfer function and decoherence from 
atmospheric seeing was performed by interleaving science 
observations with measurements of unresolved calibrators 
stars. At the end of the pipeline, the data products are 
purely interferometric as if obtained with a long-baseline 
interferometer. A full description of this experiment can be 
found in Tuthill et al. (2000b) and Monnier (1999), with 
further discussion of systematic errors in Monnier et al. 
(2004) and Monnier et al. (2007). All V 2 and closure phase 
data are available from the authors; all data products 
are stored in the FITS-based, optical interferometry data 
exchange format (OI-FITS), as described in Pauls et al. 
(2005). 
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2.3.2 Basic Results 

Before undertaking radiative transfer modeling, we provide 
the results of basic geometrical analysis of the visibility 
data. The simplest representation of the data is generally 
a circularly-symmetric Gaussian envelope, a useful model 
to give a characteristic size to the emission. Table 5 pro- 
vides the visibility intercept (Vo) and the Full- width at Half- 
maximum (FWHM) for the best fit for all datasets, including 
the reduced y 2 • Errors are generally dominated by systemat- 
ics related to the calibration procedure (i.e., seeing variation 
between source and calibrator visits) and we have used the 
relations established in Monnier et al. (2007) to quantify our 
errors. In some cases, there was evidence of two components 
to the visibility curve and we have also fitted a slightly more 
complex model of a point source plus a Gaussian envelope 
to all epochs. Table 6 contains the best fit parameters of 
the 2-component model, including the estimated fraction of 
light in the point source (/point) and the fraction of light in 
the Gaussian envelope (/Gauss)- 

In addition, we fitted each object with a 2-dimensional 
Gaussian function in order to search for signs of asymme- 
try. Objects with observed asymmetry are marked with an 
asterisk in Table 5 . Table 7 lists all the object with con- 
firmed asymmetries and we include the amount of elongation 
( FWHM maJ ° r ) an d position angle (degrees East of North) 
of the major axis. Here we have used the spread of measured 
position angles between wavelength channels and epochs to 
estimate the PA error. We will discuss further these findings 
in §4. 

3 DUST SHELL MODELING 

3.1 Introduction 

The objects in our study all have spectral energy distribu- 
tions that peak in the infrared. Indeed, these stars are sur- 
rounded by dust shells that absorb the stellar light and then 
reemit the energy in the infrared. In order to extract physical 
characteristics of these dust shells (i.e., optical depths, tem- 
peratures, etc), we must be able to compute how the dust 
will absorb, scatter, and reemit the energy from the star. We 
accomplish this with the radiative-transfer model DUSTY 
(Ivezic et al. 1999). While DUSTY is limited to calculations 
in spherical symmetry, we established in the previous section 
that most of our objects show only mild signs of global asym- 
metries; however, we caution that our results will be suspect 
for the most asymmetric of the targets listed in Table 7. 
Given a small number of input parameters, DUSTY can 
quickly compute synthetic photometry and intensity pro- 
files for dust shells. These outputs can then be compared to 
the data that we have experimentally obtained. 

3.2 Model Description 

We applied a uniform procedure for fitting all of our objects. 
Here we discuss which properties were held fixed and how 
we explored a grid of the key dust shell parameters. 

We begin with the central star. At the beginning of 
our study we used a featureless Planck blackbody spec- 
trum, however we came to realise that a blackbody spec- 
trum is a rather poor approximation for the extremely late- 


type giants in our sample due to strong molecular absorp- 
tion bands. Most notably, the HCN absorption feature of 
carbon-rich stars sits directly at the PAHcs (3.0825 fim) 
wavelength, where we have many observations. Because of 
the severe optical absorption of the dust, spectral types are 
not known for most stars in our sample and we have adopted 
an effective temperature of 2600K for all stars, which is 
as cool as we could find converged synthetic spectra. For 
the carbon stars we used a MARCS model as described in 
Loidl et al. (2001) and for the M-giants we used a PHOENIX 
NEXTGEN model as described in Hauschildt et al. (1999). 
The medium-resolution synthetic spectra from these sources 
were smoothed before input into DUSTY. Unfortunately, we 
do not have useful distance estimate to our sources - so we 
adopted a distance of 1000 pc and interstellar reddening 
of Eb-v = 0.5 for all objects. We note that the dust shells 
around the stars absorb nearly all of the energy from the cen- 
tral source, acting as a kind of calorimeter. Thus, while our 
2600K estimate for the central star temperature is crude, we 
expect the bolometric luminosity (for assumed d = 1000 pc) 
to be more accurate. However in practice our luminosity 
estimates are poor due to uncertainties in the dust shell op- 
tical depth and the fact we are not integrating the whole 
observed SED throughout the mid- and far-infrared. 

Based on the shape of the SED (and the presence 
of a silicate feature in IRAS-LRS spectra), we deter- 
mined each star to have either carbon-rich dust or silicate- 
rich dust. Based on this assignment, we chose amor- 
phous carbon (Hanner 1988) or warm amorphous silicates 
(Ossenkopf et al. 1992) respectively in the DUSTY model 
setup. Speck et al. (2008) discussed how silicates close to 
AGB stars could quickly anneal to crystalline grains but a 
full exploration of optical constants for different grain types 
was beyond the scope of this work. For the grain size distri- 
bution, we adopted the standard MRN power-law grain size 
distribution between 0.005-0.25 /am (Mathis et al. 1977); a 
later exploration of larger grain sizes did not systematically 
improve fits (also see discussion by Speck et al. 2009). An- 
other property of the dust shell we fixed is that the dust 
density follows a r~ 2 power- law, corresponding to constant 
mass- loss rate. 

Lastly, we come to the parameters of the model that 
are not fixed: the temperature of the dust shell at the inner 
boundary, Tdust, the radius of the star, R s tar, and the Ri- 
band optical depth T 2 . 2 ^m of the dust shell (as integrated 
along the line-of-sight from the observer to the star). In the 
next section, we explain our fitting procedure. 

3.3 Fitting Methodology 

We explored inner dust temperatures Tdust between 400K- 
1500K. This range explored both the high temperatures 
thought to be prohibitive of dust creation and low tem- 
peratures too cool for steady- state dust production. Note 
that when setting up a model in DUSTY, one does not 
specify the inner radius of the dust shell: this quantity is 
calculated based on the luminosity of the star and the spe- 
cific inner shell dust temperature Tdust- In terms of optical 
depth, we explored T 2 . 2 ^m between 0 and 9. This range pro- 
vided a full fitting region for our objects and values of T 2 . 2 ^ m 
much above 9 were too computationally expensive. Finally, 
Rstar was recognised to simply be a scaling factor for the 
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model outputs and could easily be optimised for every pair 
of (Tdust, t 2 . 2 ixm)- Because the DUSTY calculation was fast 
and we only had to optimise over a few parameters, we chose 
to carry out an exhaustive grid calculation over all (Tdust? 

T2.2/xm)- 

For each location in the grid we calculated the model 
SED as well as the radial intensity profiles. We calculated 
a x 2 based on both our coeval near-infrared photometry as 
well as Keck masking visibility curves. For the SED, we also 
used including V-band magnitudes in our fit with a very low 
weight to ensure that the optical depths were not too low 
(important especially when for objects without photometry 
in all three near-IR wavelength bands). When calculating 
the x 2 f° r the visibility curves, we adopted the following 
procedure. Because the y-intercept of our observed visibil- 
ity data can fluctuate ±5% due to seeing calibration, we 
normalised each visibility to 1.0 at zero baseline before fit- 
ting. Also, we weighted the visibility points so that the SED 
and the visibility data were separately given equal weight 
in the final reduced x 2 * We purposefully chose not to in- 
clude longer wavelength SED measurements, such as IRAS 
data, in our fitting. By fitting only to near-IR photometry 
and near-IR spatial data we can isolate and only probe dust 
emitted within the last few decades. This allows us to keep 
the model as simple as possible and enhances the validity of 
our assumption of constant mass loss rate (ie., p oc r -2 ). 

Once the grid calculation over inner dust temperature 
T dust and 72 . 2 ^ m was completed, the x 2 surface was used to 
estimate the best-fit parameters. The uncertainty estimates 
were produced by considering the region where the reduced 
X 2 was less than 2, a highly conservative criterion that re- 
flects the highly-correlated errors in our datasets. In cases 
where the best-fit x 2 is above 1, we scaled the x 2 results by 
the best-fitting value before estimating the parameter uncer- 
tainties. The best-fitting parameters and their uncertainties 
are compiled in Table 8. 

In addition to providing the fitting results in tabulated 
form, we also include here a series of figures which graphi- 
cally represent the new data, modeling results, and the x 2 
surface in our grid. These plots can be found in each of 
Figures 1-20. The first panel in each figure contains the 
observed near-IR photometry and best-fit model SED. The 
second panel in each figure contains the multi-wavelength 
visibility curves averaged azimuthally along with the model 
curves. Finally, the third panel shows the x 2 surface in the 
(Tdust, 72.2// m) plane. We have grouped all the epochs for 
the same object together so one can see the self-consistency 
in the derived dust shell parameters - indeed, consistent 
dust shell properties were recovered when fitting to different 
epochs, despite large changes in the central star luminosity 
due to pulsations. 

One of the most important results to take away from 
these panels that we clearly break the standard degener- 
acy between dust temperature and optical depth. This is 
because of our new spatial information - by measuring the 
sizes of the dust shell at various wavelengths we can simul- 
taneously constrain the temperature and optical depth. In 
the past, one typically had to choose an inner dust tem- 
perature based on physical arguments concerning the dust 
condensation temperatures of various dust species. Here, we 
see that the inner dust temperature can be constrained in- 


dependently from other parameters and the implications are 
discussed further in the next section. 

While the simultaneous fits to the near-IR SED and 
visibility data were generally acceptable, we found the fits 
to the shortest wavelength visibility data at H band were 
systematically worse. Since this band is most sensitive to 
scattering by dust, we explored modified dust distributions, 
especially using larger grains; we did not find systematic 
improvements to the fits by altering dust size distribution 
from MRN or by using other dust constants. 


4 DISCUSSION 

Our survey provides the first constraints on the asymmetry 
of the dust shells for such a large sample of dust-enshrouded 
AGB stars. We found that 4 out of 7 M-stars and 5 of 13 
C-stars showed evidence of dust shell asymmetries, with 
dust shell elongations between 10% and 40%. While this 
level of asymmetry may sound mild, it actually (quanti- 
tatively) compares to the level of asymmetry that would 
be expected for the most asymmetric dust shells known if 
placed at 1 kpc. For instance, we know that IRC +10216 
(Tuthill et al. 2000a) and CIT 6 (Monnier et al. 2000) have 
dramatic global asymmetries in their dust shell, detailed 
imaging made possible by virtue of their proximity. If we 
placed these targets farther away, we would not be able to 
image the detail but they would appear ~20% elongated, 
similar to the degree observed here in 45% of our sample. For 
CIT 3, we confirm the asymmetries seen by Hofmann et al. 
(2001) and note that Vinkovic et al. (2004) showed that the 
20% elongation could be explained by a bipolar outflow. 
That said, clumpy dust formation (Fleischer et al. 1992) 
might also cause stochastic variations in the inner dust shell 
geometry that could appear as short-lived elongations. Mid- 
infrared observations with long-baseline interferometers (e.g, 
ISI, VLTI-MIDI) should focus on these targets to determine 
the nature of the asymmetries. In addition, long-term moni- 
toring of these dust shells will help settle debates concerning 
when the environments of evolved stars develop large scale 
asymmetries commonly revealed in the later planetary neb- 
ula stage. For instance, a long-term asymmetry in a constant 
position angle (as judged by linear polarization or spatially 
resolved data) would be a sign of a global bipolar mass-loss 
asymmetry and not just weather. 

In order to look at dust shell properties for our full 
sample, we have plotted the inner edge dust temperature 
Tdust vs total dust shell optical depth 72 . 2 ^™, for all our 
targets. Figure 21 shows these results split into O-rich 
and C-rich dust type. For K-band optical depths below 2, 
we find the sublimation temperature of T su b (silicates) = 
1130+90 K and T su b (amorphous carbon) = 1170+60 K, both 
somewhat lower than expected from laboratory measure- 
ments (Lodders & Fegley 1999) and vastly below temper- 
atures inferred from the inner edge of YSO disks (~1800K, 
Tannirkulam et al. 2008; Benisty et al. 2010). One compo- 
nent to the observed lower dust temperature could be due 
to the fact that the central star varies in luminosity by 
about a factor of 2 during the pulsation cycle and we see 
the dust cooler than the condensation temperature during 
phases away from maximum light. 

The Tdust vs optical depth T 2 . 2 fj.rn diagram (Figure 21) 
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also shows no statistically-significant difference between O- 
rich and C-rich dust types, counter to expectation of higher 
temperatures for carbon- rich dust (Lodders V Fegley 1999). 
We recognise that our simple dust shell modeling may not 
lead to accurate estimates of the dust sublimation temper- 
ature if the inner dust formation environment radically de- 
parts from a power law density distribution, perhaps due 
to pulsations, timescale for dust formation, or multiple dust 
species. Interestingly though these concerns would likely af- 
fect C-rich and O-rich shells similarly and so the lack of a 
clear difference in sublimation temperatures between these 
dust types appears robust. 

The other important feature of Figure 21, Tdust vs op- 
tical depth T2.2/J, m 5 is the apparent temperature at the inner 
edge of the dust shell gets lower and lower with increasing 
optical depths above 2. This appears true for both C-rich 
and O-rich shells. Here we do not believe we are seeing an 
actual reduction in the dust sublimation temperature, but 
rather a change in the temperature profile in the inner dust 
formation zone due to a breakdown in the assumption of a 
spherically-symmetric r~ 2 density power law. We have am- 
ple evidence that dust formation is clumpy, as has been im- 
aged in great detail for IRC +10216 (Tuthill et al. 2000a), 
but these clumps have been shown to have a relatively weak 
affect on the temperature structure for low optical depths. 
Next we further explore how a clumpy dusty environment 
could change the temperature profile of the dust shell when 
the individual clumps become themselves optically thick to 
the stellar and even hot dust radiation field. 

Clumpy structures are seen to evolve in 2-D models 
of dust shells due to self-amplifying density perturbations 
(e.g. Woitke et al. 2000). First optically thick dust regions 
form and these regions cast shadows on the dust behind 
them. Consequently, the temperatures decrease by 100’s of 
degrees K and this allows for a higher rate of dust forma- 
tion in these shadow regions. Scattering and re-emission of 
light by the optical thick regions increases the intensity of 
radiation between them and eventually the light escapes 
through the optically thin regions in between the optically 
thick regions. Thereupon, the temperature within the opti- 
cally thin regions increases, which decreases the rate of dust 
production. These processes thus amplify the initial homo- 
geneities until large-scale clumpy structures start to form, 
such as “dust fingers” (Woitke V Niccolini 2005). Indeed, 
Woitke & Niccolini (2005) did see average dust temperatures 
to be reduced due to these opacity effects but at much weaker 
level than we see in Figure 21. Realizing that our data reveal 
a strong effect only at r’s several times larger than probed by 
Woitke V Niccolini (2005), we suggest that dust shadowing 
effects get dramatically stronger when individual clumps be- 
come optically thick to both stellar radiation as well as hot 
dust emission. A 3D radiative transfer calculation of a dusty 
dust shell could validate or disprove this explanation. 

In conclusion, our large sample of spatially resolved 
dust-enshrouded stars have led to new insights into the late 
stages of AGB star evolution. We find levels of dust shell 
elongations that point to significant asymmetries in nearly 
half of our targets. Our spatial and SED data combined 
has eliminated some model degeneracies, and we now have 
the best constraints on the actual sublimation temperatures 
for dust forming in this outflows, finding lower temperatures 
than expected from terrestrial experiments and not confirm- 


ing the large difference expected between carbon-rich and 
silicate-rich dust. Lastly, we discovered a systematic change 
in the temperature profile for inner-most dust regions when 
the dust shell optical depth rises above T 2 . 2 ^m > 2. This 
observed lowering of the central dust temperatures could be 
naturally explained as a consequence of shadowing caused 
by clumpy dust formation on spatial scales smaller than our 
angular resolution, but other possibilities should be further 
explored as well. 
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Table 1 . Properties of NIRC Camera Infrared Filters. Reference: The NIRC Manual. 


Name 

Center Wavelength 
Ao (^ m ) 

Bandpass FWHM 
AA ( fi m) 

Fractional 

Bandwidth 

Fell 

1.6471 

0.0176 

i.i% 

H 

1.6575 

0.333 

20% 

K 

2.2135 

0.427 

19% 

Kcont 

2.25965 

0.0531 

2.3% 

CH4 

2.269 

0.155 

6.8% 

PAHcs 

3.0825 

0.1007 

3.3% 


Table 2. Basic Properties of Targets 


Source 

RA (J2000) 

Dec (J2000) 

V 

J a 

H a 

K a 

Spectral 

Names 



mag 

mag 

mag 

mag 

Type 

AFGL 230 

01 33 51.21 

+62 26 53.5 

— 

16.747 

11.232 

7.097 

mV) 

AFGL 2019 

17 53 18.9 

-26 56 37 

20.2( 2 ) 

6.338 

4.035 

2.616 

M8W 

AFGL 2199 

18 35 46.48 

+05 35 46.5 

— 

8.04 

4.85 

2.701 

M( 0 ) 

AFGL 2290 

18 58 30.02 

+06 42 57.7 

— 

13.169 

8.966 

5.862 

M( 0 ) 

CIT 1 

00 06 52.94 

+43 05 00.0 

9.00^) 

3.041 

1.829 

1.115 

M9W 

CIT 3 

01 06 25.98 

+12 35 53.0 

— - 

7.45 

4.641 

2.217 

M9) 1 ) 

vl300 Aql 

20 10 27.87 

-06 16 13.6 

20^) 

6.906 

3.923 

2.059 

M« 

AFGL 1922 

17 07 58.24 

-24 44 31.1 

— 

12.244 

9.181 

6.342 

C< s > 

AFGL 1977 

17 31 54.98 

+17 45 19.7 

9 . 9 W 

10.536 

7.994 

5.607 

cm 

AFGL 2135 

18 22 34.50 

+27 06 30.2 

— 

9.043 

6.002 

3.643 

cm 

AFGL 2232 

18 41 54.39 

+17 41 08.5 

9 . 7 W 

5.742 

3.444 

1.744 

cm 

AFGL 2513 

20 09 14.22 

+31 25 44.0 

— 

8.229 

5.705 

3.69 

cm 

AFGL 2686 

20 59 08.88 

+27 26 41.7 

20*V 

9.112 

6.268 

4.075 

c e m 

AFGL 4211 

15 11 41.89 

-48 20 01.3 

— 

10.711 

7.751 

5.154 

cm 

IRAS 15148-4940 

15 18 22.05 

-49 51 04.6 

11.8W 

5.297 

3.071 

1.696 

cm 

IY Hya 

10 17 00.52 

-14 39 31.4 

14U) 

5.919 

3.666 

1.964 

cm 

LP And 

23 34 27.66 

+43 33 02.4 

— 

9.623 

6.355 

3.859 

cm 

RV Aqr 

21 05 51.68 

-00 12 40.3 

11.5W 

4.046 

2.355 

1.239 

cm 

vl899 Cyg 

21 04 14.8 

+53 21 03 

15.6W 

10.84 

8.693 

6.596 

csm 

V Cyg 

20 41 18.2702 

+48 08 28.835 

7.7W 

3.096 

1.273 

0.117 

cm 


a These magnitudes (from 2MASS) are merely representative since the targets are variable. See Table 3 for our new photometry. 
Note: Horizontal line separates oxygen-rich (top) from carbon-rich (bottom). 

References: (1) Simbad, (2) Monet (1998), (3) Buscombe (1998), (4) Egret et al. (1992), (5) Skiff (2007), (6) Volk (2002) (7) 
Garcia-Hernandez et al. (2007) 
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Table 3. Journal of observations and derived photometry 


Target 

Date(s) 

(UT) 

Filter 

Aperture 

Mask 

Magnitude 

Calibrator 

Names 

AFGL 230 

1997 Dec 

k 

FFA 

8.34±0.1 

X Cas 



PAHcs 

KL Relation* 

5.11±0.2 



2002 Jul 

k 

FFA 

8.99±0.1 

HD 9878 



PAHcs 

FFA 

5.90±0.3 

HD 9329 

AFGL 2019 

2000 Jun 

CH4 

annulus 36 

2.48±0.1 

HD 163428 



h 

annulus 36 

3.84± 0.1 

HD 156992 



PAHcs 

annulus 36 

1.52±0.1 

HD 163428 

AFGL 2199 

1998 Apr 

CH4 

annulus 36 

2.99±0.1 

HD 170137 



PAHcs 

annulus 36 

1.80±0.1 

HD 170137 

AFGL 2290 

1998 Jun 

CH4 

annulus 36 

4.72±0.1 

HD 173074 



PAHcs 

annulus 36 

2.60±0.1 

HD 173074 


1999 Apr 

CH4 

annulus 36 

5.61±0.1 

HD 173833 



k 

annulus 36 

6.19±0.32 

HD 231437 



PAHcs 

annulus 36 

3.29±0.1 

HD 173833 

CIT 1 

2000 Jun 

CH4 

annulus 36 

2.60±0.1 

A And 



h 

annulus 36 

4.18±0.25 

HD 222499 



PAHcs 

annulus 36 

1.51±0.1 

A And 

CIT 3 

1997 Dec 

Kcont 

annulus 36 

1.08±0.1 

<5 Psc 



PAHcs 

annulus 36 

-0.14±0.1 

<5 Psc 


1998 Sep 

CH4 

Golay 21 

2.45±0.1 

5 Psc 



PAHcs 

Golay 21 

1.04±0.1 

<5 Psc 

vl300 Aql 

1998 Jun 

CH4 

annulus 36 

1.39±0.1 

HD 189114 



h 

annulus 36 

3.29±0.25 

HD 192464 



PAHcs 

annulus 36 

0.60±0.1 

HD 189114 


1999 Jul 

kcont 

annulus 36 

2.02±0.1 

SAO 14382 



PAHcs 

annulus 36 

0.86±0.1 

SAO 14382 

AFGL 1922 

2000 Jun 

k 

annulus 36 

6.34±0.25 

HD 156992 



PAHcs 

KL relation* 

3.62±0.25 



2001 Jun 

k 

annulus 36 

4.90±0.1 

HD 158774 



PAHcs 

KL relation* 

2.32±0.25 


AFGL 1977 

1998 Jun 

CH4 

annulus 36 

4.19±0.1 

HD 158227 



h 

annulus 36 

7.05±0.1 

HD 158227 



PAHcs 

annulus 36 

1.84±0.1 

HD 157049 


1999 Apr 

CH4 

annulus 36 

2.77±0.1 

HD 157049 



PAHcs 

annulus 36 

0.59±0.1 

HD 157049 

AFGL 2135 

2001 Jun 

k 

annulus 36 

3.29±0.1 

HD 168366, HD 181700 



PAHcs 

annulus 36 

1.27±0.3 

HD 177716 

AFGL 2232 

1998 Jun 

CH4 

annulus 36 

2.04±0.1 

HD 158227 



h 

annulus 36 

4.12±0.1 

HD 158227 



PAHcs 

annulus 36 

0.68±0.1 

HD 157049 



CH4 

Golay 21 

2.28±0.3 

HD 168720 



PAHcs 

Golay 21 

0.94±0.3 

HD 168720 


1999 Apr 

CH4 

annulus 36 

1.06±0.1 

HD 173833 



PAHcs 

annulus 36 

-0.38±0.1 

HD 173833 

AFGL 2513 

1998 Sep 

h 

annulus 36 

6.58±0.1 

HD 196241 



CH4 

annulus 36 

4.03±0.1 

HD 200451 



PAHcs 

annulus 36 

3.16±0.3 

e Cyg 


1999 Jul 

CH4 

annulus 36 

2.90±0.1 

HD 188947 



PAHcs 

annulus 36 

1.70±0.1 

HD 188947 

AFGL 2686 

1998 Sep 

CH4 

annulus 36 

2.95±0.1 

HD 200451 



h 

annulus 36 

5.82±0.21 

HD 200451 



PAHcs 

annulus 36 

l.OliO.l 

e Cyg, A And 


1999 Jul 

CH4 

annulus 36 

5.13±0.1 

HD 188947 



PAHcs 

annulus 36 

2.92±0.1 

HD 188947 



h 

annulus 36 

8.48±0.3 

HD 198330 

AFGL 4211 

2000 Jun 

CH4 

annulus 36 

3.62±0.3 

HD 137709 



PAHcs 

annulus 36 

1.42±0.3 

HD 137709 


2001 Jun 

k 

annulus 36 

4.70±0.1 

HD 137709 



PAHcs 

KL relation* 

2.64±0.2 


IRAS 15148-4940 

2001 Jun 

CH4 

annulus 36 

1.25±0.3 

HD 137709 



k 

annulus 36 

1.30±0.3 

HD 137709 



PAHcs 

annulus 36 

1.71±0.1 

HD 136422 

IY Hya 

1999 Apr 

CH4 

annulus 36 

2.08±0.1 

HD 87262 



PAHcs 

annulus 36 

1.37±0.1 

ft Hya 

LP And 

1998 Sep 

CH4 

annulus 36 

3.89±0.1 

HD 222499, A And 



h 

annulus 36 

7.05±0.25 

HD 222499 



PAHcs 

annulus 36 

1.72±0.1 

A And 


1999 Jul 

CH4 

Golay 21 

4.01±0.1 

a Cas 



PAHcs 

Golay 21 

1.80±0.1 

a. Cas 


1999 Jan 

CH4 

Golay 21 

3.26±0.1 

a Cas 



PAHcs 

Golay 21 

1.18±0.1 

a Cas 

RV Aqr 

1999 Jul 

CH4 

Golay 21 

1.23±0.25 

SAO 143482, 3 Aqr 



PAHcs 

Golay 21 

0.56±0.25 

SAO 143482, 3 Aqr 


1998 Jun 

CH4 

Golay 21 

1.52±0.1 

HD 196321 



PAHcs 

Golay 21 

1.15±0.1 

HD 196321 

vl899 Cyg 

1998 Jun 

CH4 

annulus 36 

5.53±0.1 

HD 202897 



h 

annulus 36 

7.87±0.3 

HD 200817 



PAHcs 

annulus 36 

3.71±0.1 

HD 202897 


1999 Jul 

k 

annulus 36 

6.40±0.1 

HD 198661 



PAHcs 

KL relation* 

4.72±0.2 


V Cyg 

1998 Jun 

feii 

annulus 36 

2.59±0.1 

HD 192909 



kcont 

annulus 36 

0.53±0.1 

HD 192909 



CH4 

Golay 21 

0.50±0.1 

HD 192909 



PAHcs 

annulus 36 

0.26±0.1 

HD 192909 



PAHcs 

Golay 21 

0.19±0.1 

HD 192909 


1999 Apr 

CH4 

Golay 21 

0.15±0.1 

€ Cyg 



PAHcs 

Golay 21 

-0.25±0.1 

£ Cyg 


2001 Jun 

CH4 

annulus 36 

-0.27±0.1 

£ Cyg 



PAHcs 

annulus 36 

-0.69±0.1 

£ Cyg 


* This point was extrapolated 
Note: Horizontal line separates 


from another epoch for the same star and assigned 
oxygen-rich (top) from carbon-rich (bottom). 


error of 0.2 mag. 
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Table 4. Basic Properties of Calibrators 


Calibrator 

J 

mag 

H 

mag 

K 

mag 

PAHcs 

mag 

Reference 

HD 168720 

1.79 

0.875 

0.870 

0.794 

McWilliam Sz Lambert (1984), Cutri (2003), Neugebauer Sz Leighton (1969) 

HD 170137 

3.476 

2.737 

2.230 

2.16 

Cutri (2003), Neugebauer Sz Leighton (1969) 

e Cyg 

0.641 

0.2 

0.1 

0.011 

Neugebauer Sz Leighton (1969), Ghosh et al. (1984), Price Sz Murdock (1999) 

HD 200451 

4.101 

3.231 

2.840 

— 

Cutri (2003), Neugebauer Sz Leighton (1969) 

HD 231437 

5.027 

3.958 

3.693 

— 

Cutri (2003) 

HD 173833 

3.488 

2.647 

2.1 

2.02 

Cutri (2003), Neugebauer Sz Leighton (1969) 

HD 158227 

5.626 

4.984 

4.812 

— 

Cutri (2003) 

HD 157049 

1.975 

1.149 

.830 

.684 

Cutri (2003), Neugebauer Sz Leighton (1969), Price Sz Murdock (1999) 

HD 168366 

5.049 

4.535 

4.255 

— 

Cutri (2003) 

HD 181700 

3.938 

2.993 

2.735 

— 

Cutri (2003) 

SAO 143482 

1.665 

0.790 

0.573 

0.436 

Cutri (2003), Gullixson et al. (1983) 

HD 189114 

3.212 

2.030 

1.953 

1.908 

Cutri (2003), Gosnell et al. (1979) 

HD 137709 

2.232 

1.532 

1.331 

1.257 

Cutri (2003), extrapolation 

HD 222499 

4.641 

3.804 

3.627 

— 

Cutri (2003) 

A And 

1.970 

1.4 

1.287 

1.245 

Johnson et al. (1966), Price Sz Murdock (1999), Selby et al. (1988) 

HD 9878 

6.631 

6.730 

6.698 

— 

Cutri (2003) 

HD 9329 

4.961 

4.381 

4.341 

4.29 

Cutri (2003), extrapolation 

HD 156992 

3.901 

3.123 

2.926 

— 

Cutri (2003) 

HD 158774 

4.403 

3.451 

3.138 

— 

Cutri (2003), Kawara et al. (1983) 

HD 198611 

3.755 

2.862 

2.470 

— 

Cutri (2003), Cutri (2003), Neugebauer Sz, Leighton (1969) 

HD 202987 

3.859 

3.067 

2.82 

2.75 

Cutri (2003), Neugebauer Sz Leighton (1969) 

3 Aqr 

0.934 

-0.020 

-0.220 

-0.338 

Carter (1990) 

HD 192909 

1.190 

— 

0.180 

0.101 

Johnson et al. (1966), Neugebauer Sz Leighton (1969), Price Sz Murdock (1999) 

C c yg 

.995 

0.130 

-0.070 

-0.150 

Johnson et al. (1966), Noguchi et al. (1981) 

HD 200817 

4.174 

3.721 

3.708 

— 

Cutri (2003) 

HD 192464 

5.180 

4.176 

3.879 

— 

Cutri (2003) 

a Cas 

0.371 

-0.191 

-0.270 

-0.399 

Voelcker (1975), Alonso et al. (1994) 

li Hya 

1.216 

.506 

0.37 

0.28 

Cutri (2003), Price Sz Murdock (1999), Johnson et al. (1966) 

HD 87262 

2.974 

2.052 

1.880 

— 

Cutri (2003), Price Sz Murdock (1999), Neugebauer Sz Leighton (1969) 

HD 196321 

2.128 

1.361 

1.21 

.98496 

Cutri (2003), Price Sz Murdock (1999), Neugebauer Sz Leighton (1969) 

HD 136422 

— 

— 

0.8 

0.535 

Price (1968), Price Sz Murdock (1999), Eggen (1969) 

S Psc 

2.031 

1.198 

.890 

0.739 

Cutri (2003), Gosnell et al. (1979) 

HD 198330 

4.988 

4.159 

3.816 

— 

Cutri (2003) 

HD 188947 

1.934 

1.438 

1.621 

1.561 

Noguchi et al. (1981), Elias et al. (1982), Glass (1975) 

X Cas 

3.019 

2.481 

2.311 


Cutri (2003), Neugebauer Sz Leighton (1969) 

HD 163428 


— 

1.6 

1 . Hi 1 

White Sz Wing (1978), Humphreys Sz Ney (1974) 

HD 196241 

1.19 

3.620 

3.090 

— 

Morel Sz Magnenat (1978), Cutri (2003) 
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Table 5. Results from Circular ly-Symmetric Gaussian Models 


Target 

Date(s) 

Filter 

Aperture 

v 0 

FWHM 

x 2 /dof 





(±0.05) 

(mas) 


AFGL 230 

1997 Dec 

k 

FFA 

0.71 

32±3 

0.23 


2002 Jul 

k 

FFA 

0.54 

34±3 

0.34 



PAHcs 

FFA 

0.74 

33±2 

0.05 

AFGL 2019 

2000 Jun 

CH4 

annulus 36 

0.96 

10 + 6 
±u -10 

0.31 



h 

annulus 36 

0.90 

9±4 

0.65 



PAHcs 

annulus 36 

0.95 

21±3 

0.27 

AFGL 2199 

1998 Apr 

CH4 

annulus 36 

0.92 

14±6 

0.23 



PAHcs 

annulus 36 

1.00 

22±3 

0.45 

AFGL 2290* 

1998 Jun 

CH4 

annulus 36 

0.76 

22±4 

0.33 



PAHcs 

annulus 36 

0.84 

27±3 

0.69 


1999 Apr 

CH4 

annulus 36 

0.72 

34±3 

0.39 



k 

annulus 36 

0.75 

32±3 

0.51 



PAHcs 

annulus 36 

0.83 

36±2 

0.36 

Cit 1* 

2000 Jun 

CH4 

annulus 36 

0.92 

15±5 

0.35 



h 

annulus 36 

0.93 

14±3 

0.37 



PAHcs 

annulus 36 

0.94 

20±4 

0.46 

Cit 3* 

1997 Dec 

kcont 

annulus 36 

0.89 

20±5 

0.44 



PAHcs 

annulus 36 

0.89 

37±2 

0.21 


1998 Sep 

CH4 

Golay 21 

0.89 

21±4 

0.35 



PAHcs 

Golay 21 

0.90 

29±2 

0.25 

vl300 Aql* 

1998 Jun 

CH4 

annulus 36 

0.83 

14±6 

0.43 



h 

annulus 36 

0.81 

14±3 

0.41 



PAHcs 

annulus 36 

0.84 

23±3 

0.50 


1999 Jul 

kcont 

annulus 36 

0.87 

18±5 

0.39 



PAHcs 

annulus 36 

0.90 

21±3 

0.52 

AFGL 1922 

2000 Jun 

k 

annulus 36 

0.76 

24±4 

0.88 


2001 Jun 

k 

annulus 36 

0.83 

29±4 

0.76 



PAHcs 

annulus 36 

0.95 

58±2 

0.43 

AFGL 1977* 

1998 Jun 

CH4 

annulus 36 

0.78 

24±4 

0.26 



h 

annulus 36 

0.76 

17±3 

0.68 



PAHcs 

annulus 36 

0.94 

34±2 

0.41 


1999 Apr 

CH4 

annulus 36 

0.96 

29±4 

0.25 



PAHcs 

annulus 36 

0.89 

52±2 

0.26 

AFGL 2135 

2001 Jun 

k 

annulus 36 

0.66 

17±5 

0.49 


2001 Jun 

PAHcs 

annulus 36 

0.50 

34±2 

0.13 

AFGL 2232* 

1998 Jun 

CH4 

annulus 36 

0.83 

18±5 

1.32 



h 

annulus 36 

0.81 

14±3 

0.50 



PAHcs 

annulus 36 

0.90 

33±2 

0.56 



CH4 

Golay 21 

0.91 

20±5 

0.17 



PAHcs 

Golay 21 

0.90 

34±2 

0.14 


1999 Apr 

CH4 

annulus 36 

0.69 

44±3 

0.19 



PAHcs 

annulus 36 

0.86 

42±2 

0.12 

AFGL 2513* 

1998 Sep 

h 

annulus 36 

1.00 

ill 

1.15 



CH4 

annulus 36 

0.94 

0 

1 + 
M C5 

o 

0.18 



PAHcs 

annulus 36 

1.00 

16±4 

0.70 


1999 Jul 

CH4 

annulus 36 

1.00 

H±6 

0.32 



PAHcs 

annulus 36 

0.96 

24±3 

0.36 

AFGL 2686 

1998 Sep 

CH4 

annulus 36 

0.89 

29±4 

0.44 



h 

annulus 36 

0.89 

26±2 

0.47 



PAHcs 

annulus 36 

0.89 

35±2 

0.36 


1999 Jul 

CH4 

annulus 36 

0.92 

26±4 

0.68 



PAHcs 

annulus 36 

0.91 

33±2 

0.44 



h 

annulus 36 

0.77 

28±2 

0.87 

AFGL 4211 

2000 Jun 

CH4 

annulus 36 

0.78 

31±3 

0.48 



PAHcs 

annulus 36 

0.82 

70±3 

0.10 


2001 Jun 

k 

annulus 36 

0.62 

20±5 

0.63 

IRAS 15148-4940 

2001 Jun 

CH4 

annulus 36 

0.77 

13±7 

0.41 



k 

annulus 36 

0.82 

13±7 

0.59 



PAHcs 

annulus 36 

0.86 

25±3 

0.39 

IY Hya 

1999 Apr 

CH4 

annulus 36 

0.88 

14±6 

0.32 



PAHcs 

annulus 36 

0.94 

33±2 

0.28 

LP And* 

1998 Sep 

CH4 

annulus 36 

0.83 

25±4 

0.52 



h 

annulus 36 

0.68 

20±3 

0.56 



PAHcs 

annulus 36 

0.86 

47±2 

0.49 


1999 Jul 

CH4 

Golay 21 

0.89 

24±4 

0.99 



PAHcs 

Golay 21 

0.79 

48±2 

0.50 


1999 Jan 

CH4 

Golay 21 

0.70 

25±4 

2.41 



PAHcs 

Golay 21 

0.66 

35±2 

0.45 

RV Aqr* 

1999 Jul 

CH4 

Golay 21 

1.00 

8±8 

0.16 



PAHcs 

Golay 21 

0.96 

26±3 

0.21 


1998 Jun 

CH4 

Golay 21 

0.98 

12±8 

0.13 



PAHcs 

Golay 21 

1.00 

27±3 

0.36 

vl899 Cyg 

1998 Jun 

CH4 

annulus 36 

0.88 

18±5 

0.35 



h 

annulus 36 

0.86 

16±3 

0.39 



PAHcs 

annulus 36 

0.92 

22±3 

0.32 


1999 Jul 

k 

annulus 36 

0.93 

15±5 

0.62 

V Cyg 

1998 Jun 

feii 

annulus 36 

0.87 

14±3 

0.92 



kcont 

annulus 36 

0.96 

16±5 

1.15 



PAHcs 

annulus 36 

0.90 

34±2 

0.51 



CH4 

Golay 21 

1.00 

18±5 

0.35 



PAHcs 

Golay 21 

0.92 

38±2 

0.10 


1999 Apr 

CH4 

Golay 21 

0.93 

17±5 

0.15 



PAHcs 

Golay 21 

0.86 

38±2 

0.12 


2001 Jun 

CH4 

annulus 36 

0.82 

19±5 

0.26 



PAHcs 

annulus 36 

0.83 

42±2 

0.15 


* target is asymmetric, see Table 7 for further details 

Note: Horizontal line separates oxygen-rich (top) from carbon-rich (bottom). 
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Table 6. Results from Central Point plus Circular ly-Symmetric Gaussian Models 


Target 

Date(s) 

Filter 

Aperture 

fpoint 

(±.05) 

^Gauss 

(±0.05) 

FWHM 

(mas) 

x 2 /dof 

AFGL 230 

1997 Dec 

k 

FFA 

0.24 

0.52 

47±7 

0.21 


2002 Jul 

k 

FFA 

0.30 

0.42 

98±7 

0.15 



PAHcs 

FFA 

0.50 

0.32 

86±8 

0.42 

AFGL 2019 

2000 Jun 

CH4 

annulus 36 

0.86 

0.14 

51 + 30 
51 -45 

0.28 



h 

annulus 36 

0.00 

0.83 


0.83 



PAHcs 

annulus 36 

0.45 

0.51 

31±5 

0.27 

AFGL 2199 

1998 Apr 

CH4 

annulus 36 

0.38 

0.54 

19±9 

0.22 



PAHcs 

annulus 36 

0.36 

0.68 

30±4 

0.43 

AFGL 2290* 

1998 Jun 

CH4 

annulus 36 

0.44 

0.38 

46± 12 

0.26 



PAHcs 

annulus 36 

0.55 

0.38 

68±7 

0.61 


1999 Apr 

CH4 

annulus 36 

0.31 

0.51 

66±9 

0.21 



k 

annulus 36 

0.34 

0.56 

66±9 

0.31 



PAHcs 

annulus 36 

0.25 

0.60 

49±3 

0.35 

CiT l* 

2000 Jun 

CH4 

annulus 36 

0.00 

0.92 

14±6 

0.35 



h 

annulus 36 

0.47 

0.49 

23±6 

0.36 



PAHcs 

annulus 36 

0.00 

0.94 

20±4 

0.46 

CIT 3 * 

1997 Dec 

kcont 

annulus 36 

0.58 

0.50 

53± 13 

0.23 



PAHcs 

annulus 36 

0.36 

0.62 

60±4 

0.02 


1998 Sep 

CH4 

Golay 21 

0.50 

0.44 

40± 10 

0.03 



PAHcs 

Golay 21 

0.46 

0.47 

50±5 

0.20 

vl300 Aql* 

1998 Jun 

CH4 

annulus 36 

0.64 

0.23 

43±20 

0.40 



h 

annulus 36 

0.45 

0.38 

25±7 

0.39 



PAHcs 

annulus 36 

0.64 

0.34 

80± 10 

0.30 


1999 Jul 

kcont 

annulus 36 

0.00 

0.87 

18±5 

0.39 



PAHcs 

annulus 36 

0.17 

0.73 

24±4 

0.52 

AFGL 1922 

2000 Jun 

k 

annulus 36 

0.42 

0.39 

47±11 

0.83 


2001 Jun 

k 

annulus 36 

0.43 

0.49 

57±12 

0.62 



PAHcs 

annulus 36 

0.51 

0.47 

105±6 

0.41 

AFGL 1977* 

1998 Jun 

CH4 

annulus 36 

0.36 

0.45 

41±9 

0.22 



h 

annulus 36 

0.50 

0.33 

43± 8 

4 -13 

0.59 



PAHcs 

annulus 36 

0.42 

0.58 

56±4 

0.36 


1999 Apr 

CH4 

annulus 36 

0.34 

0.69 

43±7 

0.17 



PAHcs 

annulus 36 

0.25 

0.74 

45±3 

0.17 

AFGL 2135 

2001 Jun 

k 

annulus 36 

0.00 

0.66 

17±5 

0.49 


2001 Jun 

PAHcs 

annulus 36 

0.28 

0.27 

78±6 

0.10 

AFGL 2232* 

1998 Jun 

CH4 

annulus 36 

0.56 

0.33 

44± 14 

0.47 



h 

annulus 36 

0.00 

0.81 

14±3 

1.32 



PAHcs 

annulus 36 

0.50 

0.49 

66±6 

0.49 



CH4 

Golay 21 

0.52 

0.45 

39± 10 

0.10 



PAHcs 

Golay 21 

0.35 

0.60 

51±4 

0.09 


1999 Apr 

CH4 

annulus 36 

0.23 

0.56 

72±7 

0.08 



PAHcs 

annulus 36 

0.26 

0.64 

58±3 

0.07 

AFGL 2513* 

1998 Sep 

h 

annulus 36 

0.03 

1.00 

it? 

1.09 



CH4 

annulus 36 

0.81 

0.14 

39±39 

0.17 



PAHcs 

annulus 36 

0.03 

1.00 

19±4 

0.68 


1999 Jul 

CH4 

annulus 36 

0.00 

0.92 

1 ±13 
— 1 

0.65 



PAHcs 

annulus 36 

0.67 

0.37 

65±9 

0.24 

AFGL 2686 

1998 Sep 

CH4 

annulus 36 

0.30 

0.63 

42±7 

0.39 



h 

annulus 36 

0.24 

0.69 

35±4 

0.43 



PAHcs 

annulus 36 

0.39 

0.56 

58±4 

0.29 


1999 Jul 

CH4 

annulus 36 

0.21 

0.73 

32±5 

0.67 



PAHcs 

annulus 36 

0.29 

0.64 

45±3 

0.42 



h 

annulus 36 

0.63 

0.30 

138 + g 5 

0.77 

AFGL 4211 

2000 Jun 

CH4 

annulus 36 

0.36 

0.54 

59± 10 

0.27 



PAHcs 

annulus 36 

0.16 

0.70 

89±4 

0.06 


2001 Jun 

k 

annulus 36 

0.47 

0.28 

86 t? 9 

0.45 

IRAS 15148-4940 

2001 Jun 

CH4 

annulus 36 

0.62 

0.18 

44±23 

0.39 



k 

annulus 36 

0.52 

0.31 

25±13 

0.59 



PAHcs 

annulus 36 

0.52 

0.38 

49±7 

0.37 

IY Hya 

1999 Apr 

CH4 

annulus 36 

0.00 

0.88 

14±6 

0.32 



PAHcs 

annulus 36 

0.00 

0.94 

32±2 

0.28 

LP And* 

1998 Sep 

CH4 

annulus 36 

0.42 

0.49 

49± 10 

0.41 



h 

annulus 36 

0.37 

0.34 

39±9 

0.52 



PAHcs 

annulus 36 

0.30 

0.65 

74±4 

0.35 


1999 Jul 

CH4 

Golay 21 

0.47 

0.49 

47±11 

0.91 



PAHcs 

Golay 21 

0.35 

0.51 

85±5 

0.42 


1999 Jan 

CH4 

Golay 21 

0.47 

0.49 

47 ±11 

0.91 



PAHcs 

Golay 21 

0.35 

0.51 

85±5 

0.42 

RV Aqr* 

1999 Jul 

CH4 

Golay 21 

0.00 

0.96 

-i ±13 
— 1 

0.30 



PAHcs 

Golay 21 

0.33 

0.64 

34±4 

0.21 


1998 Jun 

CH4 

Golay 21 

0.74 

0.25 

30±20 

0.12 



PAHcs 

Golay 21 

0.43 

0.62 

42±4 

0.33 

vl899 Cyg 

1998 Jun 

CH4 

annulus 36 

0.32 

0.57 

24±7 

0.35 



h 

annulus 36 

0.72 

0.23 

75± 15 

0.30 



PAHcs 

annulus 36 

0.65 

0.23 

57±13 

0.31 


1999 Jul 

k 

annulus 36 

0.76 

0.23 

52±24 

0.59 

V Cyg 

1998 Jun 

feii 

annulus 36 

0.31 

0.56 

18±5 

0.91 



kcont 

annulus 36 

0.51 

0.47 

26±9 

1.14 



PAHcs 

annulus 36 

0.00 

0.90 

34±2 

0.51 



CH4 

Golay 21 

0.58 

0.47 

35± 10 

0.30 



PAHcs 

Golay 21 

0.28 

0.67 

52±3 

0.06 


1999 Apr 

CH4 

Golay 21 

0.52 

0.43 

30±9 

0.14 



PAHcs 

Golay 21 

0.31 

0.60 

57±4 

0.07 


2001 Jun 

CH4 

annulus 36 

0.50 

0.36 

40± 1 1 

0.21 



PAHcs 

annulus 36 

0.28 

0.61 

63±4 

0.09 


* target is asymmetric, see Table 7 for further details 

Note: Horizontal line separates oxygen-rich (top) from carbon-rich (bottom). 
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Table 7. Results from 2-dimensional Gaussian Models 


Target 

Date(s) 

Filter 

FWHM ma j or 

FWHM minor 

<PA> 

AFGL 2290 

1998 Jun 

CH4 

1.23 

58±20 



PAHcs 

1.24 



1999 Apr 

CH4 

1.24 




k 

1.24 




PAHcs 

1.06 


CIT 1 

2000 Jun 

CH4 

1.13 

133±3 



h 

1.14 




PAHcs 

1.11 


CIT 3 

1997 Dec 

Kcont 

1.19 

151±9 



PAHcs 

1.06 



1998 Sep 

CH4 

1.28 




PAHcs 

1.04 


vl300 Aql 

1998 Jun 

CH4 

1.34 

108A13 



h 

1.14 




PAHcs 

1.19 



1999 Jul 

kcont 

1.34 




PAHcs 

1.31 


AFGL 1977 

1998 Jun 

CH4 

1.11 

71±18 



h 

1.31 




PAHcs 

1.06 



1999 Apr 

CH4 

1.10 




PAHcs 

1.21 


AFGL 2232 

1998 Jun 

CH4 

1.37 

94±10 



h 

1.83 




PAHcs 

1.19 




CH4 

1.10 




PAHcs 

1.08 



1999 Apr 

CH4 

1.22 




PAHcs 

1.05 


AFGL 2513 

1998 Sep 

h 

unresolved 




CH4 

1.5 

61±8 



PAHcs 

1.39 



1999 Jul 

CH4 

1.38 




PAHcs 

1.2 


LP And 

1998 Sep 

CH4 

1.46 

108±6 



h 

2.03 




PAHcs 

1.39 



1999 Jul 

CH4 

1.64 




PAHcs 

1.36 



1999 Jan 

CH4 

1.82 




PAHcs 

1.20 


RV Aqr 

1999 Jul 

CH4 

1.49 

122A24 



PAHcs 

1.23 



1998 Jun 

CH4 

1.35 




PAHcs 

1.22 



* Sources missing from this list were found to have circular ly-symmetric dust shells (within errors). Note: PA is the mean position 
angle of the major axis (degrees East of North) for all filters and epochs. 
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Table 8. Results from DUSTY Radiative Transfer Model 


Target 

Date(s) 

Tdust 

7’2.2/Li m 

R* 

^(1 k P c ) 

x 2 /dof 



(K) 


(mas) 

(10 3 L 0 ) 


AFGL 230 

1997 Dec 

8001“ 

4 q+0.9 
4 - y -0.7 

i e:+0.5 
4 - D -0.3 

4 c ;+ 3 - 1 
+ °-l .6 

0.26 


2002 Jul 

540l??o 

7 4 + 1 - 6 
' - 4 - 1.2 

4 1+4-5 
4 — 2.1 

qi +108 
01-24 

3.86 

AFGL 2019 

2000 Jun 

llSOtsso 

n Q9 + 0 - 23 
u - y - 0.12 

q c+0.7 

241“ 

0.54 

AFGL 2199 

1998 Apr 

iiaotiS 

1 6 +1 ' 2 
i -°-0.7 

q q+ 2.0 
0 , °— 0.6 

2 ll“ 

0.06 

AFGL 2290 

1998 Jun 

850lgQ° 

q pr+0.5 
6 -°-0.5 

q 7 + 0.8 

‘ -0.7 

261 “ 

0.33 


1999 Apr 

800111° 

4-6l° 0 : 7 5 

3 9 + 1 - 7 
°- y - 0.8 

29lH 

2.63 

CIT 1 

2000 Jun 

1190 +310 

iiyu _ 230 

1 O+0.5 
1 - z - 0.2 

q cl+ 0.8 
d-O-Q 4 

24l“ 

0.60 

CIT 3 

1997 Dec 

1HO + 230 

1 4 + 0 - 7 
i - 4 -0.5 

7 q + 1.5 
‘ — 0.6 

H 611 ? 

0.34 


1998 Sep 

102011 “ 

1 Q + 0 - 7 
4,y — 0.5 

K n+ 0.8 
°- U - 0.6 

48lH 

0.29 

vl300 Aql 

1998 Jun 

1080+1“ 

n Q9 + 0 - 46 
u - yz - 0.12 

5 8+ 1 - 0 
D - 5 - 0.6 

64li2 

0.50 


1999 Jul 

116011“* 

l.60ll;l* 

5 l + 2 -0 
°- 4 - 0.8 

49l“ 

0.11 

AFGL 1922 

2000 Jun 

850l“° 

K q+0.7 
°.9 _ 0 7 

4 1+ 1 - 1 
4 - D - 0.8 

qy +21 

1.44 


2001 Jun 

85011“ 

o q+ 0.5 
°- y -0.5 

5 2+ 1 - 2 
°' Z -1.5 

sill l 

0.39 

AFGL 1977 

1998 Jun 

9101“ 

9 e+0-1 
^•°- 0.2 

4 0 +0 ' 2 
4 - U - 0.6 

3ll| 

1.69 


1999 Apr 

9901“ 

9 pl+0.5 

Z.O _ 0 2 

6 o + 0 - 8 

O.U-0.4 

68 lg 9 

0.65 

AFGL 2135 

2001 Jun 

740 +370 

200 

q 0 + 4.2 

Q 3 +65.7 
y - 4 -4.6 

i £?7+10500 
iD ‘ -123 

5.04 

AFGL 2232 

1998 Jun 

moll“ 

i-elS:i 

c i +0.5 
°- 4 -0.5 

4811S 

0.32 


1999 Apr 

130011“ 

2 . 811 + 

9.4111 

165llf 

3.29 

AFGL 2513 

1998 Sep 

1500ll so 

1 Q + 0 - 2 
1 - y - 0.2 

i 7+0.9 
i -‘ - 0.2 

5.3lS;| 

0.57 


1999 Jul 

1110 +“° 

iiiu_ 200 

-| 9+0.9 
4-Z — 0.5 

q 1 + 1-0 
°- ± -0.4 

18ll 4 

0.30 

AFGL 2686 

1998 Sep 

moilfo 

9 q+0-5 
Z - 8 -0.5 

5 3 + 1 - 1 
° - ' 4 — 0.8 

53li5 

1.85 


1999 Jul 

8201“ 

q 9 + 0.2 
- 0.2 

q q+0.4 

o,± -0.3 

19l 4 4 

1.02 

AFGL 4211 

2000 Jun 

8801“ 

q 7 + 0.7 
' -0.5 

c 7+0.3 
D - ' -0.3 

8511 

0.73 


2001 Jun 

85011“ 

4 2+0 - 9 
^• Z -0.7 

6 1+ 3 ’ 7 
°- ± - 2.0 

71+113 
1 ± -39 

3.20 

IRAS 15148-4940 

2001 Jun 

Q40 +340 * 

y4U -170 

0 9q + 0 - 46 
U - ZO - 0.12 

4 6 +0 ' 2 

^•°-0.9 

0 

1 + 

2.47 

IY Hya 

1999 Apr 

96011“ 

n 46 +0-46 
U -^°- 0.12 

4 2 +0 ' 1 
^• z - 0.2 

3311 

0.25 

LP And 

1998 Sep 

8801“ 

q n+0-5 

9.U-0 2 

4.8+1 

43l“ 

1.49 


1999 Jul 

8201“ 

q n+0-2 
,4 - U -0.5 

4.8+1 

44ll 2 2 

1.00 


1999 Jan 

8801“ 

q 9 + 0.5 
°- Z - 0.2 

c 7+ 1 -0 
D - ' -0.7 

85lH 

1.67 

RV Aqr 

1999 Jul 

1500l! 8O 

n 46 +0-23 
U -^°-0.23 

q q+0.7 
°- 4 -0.4 

55l“ 

1.09 


1998 Jun 

H90+ 310 

±±yu - 150 

n oq+0.46 
U - ZO -0.12 

5 1+ 0 - 1 
°- i -0.7 

48+ 2 

^°-13 

0.23 

vl899 Cyg 

1998 Jun 

740 +60 

‘^ U -60 

9 q+0.5 

z.a_ o 2 

2 0 +0 - 5 
Z>U -0.3 

7 ix+3.8 
' ' D -1.8 

0.53 


1999 Jul 

600 +34 ° 

ouu -200 

9 c + 2.5 
Z.O — 1 2 

1 s+ 4 - 4 

6 1+ 66 
D - 4 -4.8 

0.12 

V Cyg 

1998 Jun 

127011“* 

0 69 +0 ' 46 * 

U.OU-o 23 

6-6±8:I 

83lg 9 

3.02 


1999 Apr 

11601“° 

n 9q+0.23 
U - ZO -0.12 

9-61°;! 

1741“ 

0.21 


2001 Jun 

127011“ 

n 46 +0-23 
U -^°-0.23 

10 . 61+1 

2121 “ 

0.24 


* This star has two regions which meet our 1-cr criteria for a best fit. The particular values shown were chosen for consistency, see the 
appropriate figure for more details. 

Note: Horizontal line separates oxygen-rich (top) from carbon-rich (bottom). 
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Figure 1 . Best fit plots for AFGL230. The first row are figures for 
the epoch Dec97 and the second row is for Jul02. The first panel 
in each row shows a fit to the SED with our new photometry in- 
cluded with errors (2MASS points are plotted as squares in each 
frame for reference). The dashed line represents the contribution 
from the star, the dotted line represents dust contribution, the 
dash-dotted line represents the contribution from scattered light, 
and the solid line is the total flux. The second panel shows our 
DUSTY fits to the visibility data for each wavelength of obser- 
vations. The third panel shows the x 2 / DOF surface, with bright 
areas showing the best-fitting region. The black contour denotes 
the l-{cr} error . 


Figure 4. Best fit plots for AFGL2290. See Fig.l caption. 
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Figure 5. Best fit plots for CIT 1. See Fig.l caption. 
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Figure 2. Best fit plots for AFGL2019. See Fig.l caption. 



Figure 3. 


Best fit plots for AFGL2199. See Fig.l caption. 




400 600 800 100012001400 
T_Dust (K) 


Figure 6. Best fit plots for CIT 3. See Fig.l caption. 
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Figure 7. Best fit plots for vl300 Aql. See Fig.l caption. For 
the 1999 Jul. epoch we chose the lower-right region as the best 
fit region because it is consistent with the best fit region for the 
1998 Jun. epoch. 




0.2 


0.0 


0 2 4 6 8 10 12 

Baseline (meters) 




400 600 800 100012001400 
T_Dust (K) 


Figure 8. Best fit plots for AFGL1922. See Fig.l caption. 
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Figure 9. 


Best fit plots for AFGL1977. See Fig.l caption. 


1Q-1U 

1 

E „ 

| 10- 11 

AFGL2135 ' "/s''' 
2001 Jun J- 

A 

1.0 

0.8 

2- 

0 ^0 2.21 jim 

A — A 3.08 jim 

8 

I 6 

AFGL2135 
2001 Jun 

|Y\ 


i 

1 0.6 
> 


« 4 


Q 

C 10' 12 


0.4 

0.2 


2 




0.0 


0 

Best x 2 5.0 


1 2345 02468 10 12 400 600 800 100012001400 

Wavelength (urn) Baseline (meters) T_Dust (K) 


Figure 10. Best fit plots for AFGL2135. See Fig.l caption. 
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See Fig.l caption. 




Figure 12. Best fit plots for AFGL2513. See Fig.l caption. 



Figure 13. 


Best fit plots for AFGL2686. See Fig.l caption. 


© 0000 RAS, MNRAS 000 , 000-000 















16 


Blasius et al. 




Figure 14. Best fit plots for AFGL4211. See Fig.l caption. 



Figure 15. Best fit plots for IRAS15148-4940. See Fig.l caption. 
We chose the lower-right region as the best fit region because in 
all other cases of multiple good fitting regions the one at low tau 
and high dust temperature was the consistent region. 
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Figure 16. 


Best fit plots for IY Hya. See Fig.l caption. 
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Figure IT. Best fit plots for LP And. See Fig.l caption. 
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Figure 18. Best fit plots for RV Aqr. See Fig.l caption. 




Figure 19. 


Best fit plots for vl899 Cyg. See Fig.l caption. 
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Figure 20. Best fit plots for V Cyg. See Fig.l caption. For 
the 1998 Jun. epoch the best fitting region was chosen to be the 
lower-right region because it is consistent with the other epochs. 



Figure 21. A plot of best fit T2.2^m versus the temperature 
at the inner edge of dust shell T^ust- Open symbols are used for 
oxygen-rich dust shells and closed symbols are used for carbon- 
rich dust shells. 
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